Isospin Blockade in Transport through Vertical Double Quantum Dots 
Bernhard Wunsch a ' 1 , David Jacob b and Daniela Pfannkuche a 

a 7. Institute of Theoretical Physics, University of Hamburg, Jungiusstr. 9, D-20355 Hamburg, Germany 
b Departamento de Fisica Aplicada, Universidad de Alicante, San Vicente del Raspeig, 03690 Alicante, Spain 



Abstract 

We study the spectrum and the transport properties of two identical, vertically coupled quantum dots in a perpendicular 
magnetic field. We find correlation-induced energy crossings in a magnetic field sweep between states differing only in 
the vertical degree of freedom. Considering the influence of a slight asymmetry between the dots caused by the applied 
source-drain voltage in vertical transport experiments these crossings convert to anticrossings accompanied by the build-up 
of charge polarization which is tunable by the perpendicular magnetic field. The polarization strongly affects the vertical 
transport through the double quantum dot and is manifest in an isospin blockade and the appearance of negative differential 
conductances in the magnetic field range where the charge localization occurs. 
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Double quantum dots have recently attracted much 
interest as possible implementations of qubits and as 
model systems to study the few-particle physics of 
molecular binding [1,2,3]. In vertically coupled dou- 
ble quantum dots (DQD) the molecular binding for a 
single electron is determined by the vertical degree of 
freedom, while the lateral degree of freedom reflects 
the physics of the single dots or quasiatoms building 
up the molecule. In this paper we describe level cross- 
ings between few-particle states differing only in parity 
along the vertical direction. These crossings are based 
on Coulomb correlations between the electrons and 
may lead to strong charge polarization, which is tun- 
able by an external perpendicular magnetic field [4]. 
Thus, these crossings mark changes of the molecu- 
lar binding which originate from many-body effects. 
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Besides its importance for quantum dot physics, the 
effect illuminates the analogy between real spin and 
isospin, the latter of which describes the vertical de- 
gree of freedom [1]. This analogy manifests itself in an 
isospin blockade of transport through the double dot 
system [5, 6, 7]. Furthermore, the localization provides 
an effective two-level system which may be useful as a 
field-tunable charge qubit. 

We model the DQD in vertical direction by two par- 
allel layers separated by a distance d. In lateral di- 
rection the electrons are confined by a rotationally- 
symmetric parabolic potential of strength Huiq. There- 
fore, the vertical degree of freedom is reduced to an ad- 
ditional spin- like degree of freedom, the isospin [1]. In 
analogy to the real spin one can define a spin algebra 
for the isospin, where the z-component I z specifies the 
vertical degree of freedom. An electron with I z — +1/2 
(I z = —1/2) is situated in the upper (lower) dot. We 
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perform an exact diagonalization of the few-particle 
Hamiltonian including tunneling between the dots as 
well as the Zeeman term and the full Coulomb inter- 
action. Thus, correlations between electrons are fully 
taken into account. The strength of tunneling t is speci- 
fied by the energy splitting between the symmetric and 
antisymmetric single particle states in z-direction, t = 
— Asas/2. Due to the symmetries of the system the 
z-component of the angular momentum, M, the mag- 
nitude of the total spin, S, and the z-component of the 
spin, Sz are conserved. If both dots are identical and no 
asymmetry between the layers is present, an additional 
symmetry leads to the conservation of the isospin par- 
ity, P = 2 N ° ■ I^P ® ■ ■ .®li N °^ specified by the quantum 
number P = ±1. The isospin parity flips all isospins, 
e.g. it moves all electron orbitals from the upper dot 
to the lower dot and vice versa. Increasing the verti- 
cal magnetic field effectively leads to a stronger lateral 
confinement of the electrons. Therefore, the Coulomb 
energy increases with the magnetic field which causes 
ground state (GS) crossings between states differing in 
angular momentum and/or spin [1]. 

In this paper we discuss energy crossings in a mag- 
netic field sweep between states differing only in parity, 
e.g. in the vertical degree of freedom [4]. These cross- 
ings originate from the different scaling with magnetic 
field of the intradot Coulomb interaction that accounts 
for the interaction between electrons on the same dot 
and interdot Coulomb interaction acting between elec- 
trons on different dots. The intradot interaction in- 
creases faster with increasing magnetic field than the 
interdot interaction, that is limited to 1/d [1]. Hence, 
the Coulomb correlations in the eigenstates of the DQD 
are magnetic field dependent. 

In the following we study parity crossings in a DQD 
containing three electrons. For small tunneling t <C 
ftwo the energy splitting between the parity eigenstates 
\P — ±1) within the same set of quantum numbers 
(M, S, S z ) is due to their different tunneling energies, 
as the occupation probabilities of symmetric and an- 
tisymmetric orbitals depends on parity. However due 
to the magnetic-field-dependent Coulomb correlations 
the favored parity may change with magnetic field. In 
the example discussed here P = — 1 is preferred for 
magnetic fields B < 7.75 T whereas P = 1 is preferred 
for strong magnetic fields, B > 7.75 T. 

Parity crossings between states with the same quan- 
tum numbers M, S and S z can have drastic effects. 




Fig. 1. Angular momentum M, total spin S and ex- 
pectation value of the z-component of the isospin (I z ) 
for the three-electron GS. The peak in (I z ) illustrates 
the charge localization that corresponds to a parity anti- 
crossing (see text). Parameters [11]: Tunneling strength: 
t = —0.059 meV, confinement: Hujq = 2.96 meV, layer sep- 
aration: d = 19.6 nm, asymmetry: V z = 5.9 * 10 — 4 meV. 

Fig. 1 illustrates the characteristics of the three- 
electron ground state (GS) in a slightly asymmetric 
DQD as function of magnetic field. An asymmetry 
between the dots is modeled by a small voltage drop 
between them. We account for it by adding the term 
Vz = V z ■ I z to the Hamiltonian, where I z = J2i=i & 
is the z-component of the total isospin and V z is as- 
sumed to be much smaller than all other parameters 
V z <C t,hioo- At around B — 7.75 T there is a sharp 
peak in the average isospin z-component, correspond- 
ing to a strong charge polarization. Right at the peak, 
two electrons occupy the lower dot and only one elec- 
tron occupies the upper one. The charge localization 
is a consequence of a parity crossing in the symmet- 
ric DQD, that turns into an anticrossing due to the 
slight asymmetry between the dots. The width of the 
peak and the energy splitting at the anticrossing are 
proportional to the strength of asymmetry [4]. For the 
calculation presented in Fig. 1 we chose the param- 
eters such, that the anticrossing occurs in the GS of 
a DQD containing three electrons. But the effect of 
charge localization is neither bound to this particular 
choice of parameters nor to three electron states. We 
found similar effects also for 5 electrons and/or in 
other subsets of quantum numbers M, S and S z . 

In the following we show that the charge localization 
drastically affects the transport through the DQD. We 
assume a vertical transport setup where the current 
and the differential conductance through the DQD are 
measured as response to a source-drain voltage V s d and 
a gate voltage Vg [3,8,10]. The source-drain voltage 
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determines the width of the transport window where 
electrons from the source reservoir can tunnel through 
the DQD to the drain reservoir, whereas the gate volt- 
age changes the energies of the DQD proportionally to 
VgN where N denotes the number of electrons in the 
DQD. The energy needed for a transition between an 
^-electron state i with energy E(N, i) and an N + 1- 
electron state j with energy E(N + l,j) is therefore 
given by fi(N + N, i) = E(N + - E(N, i) - 
eaVG, where a denotes a constant proportionality fac- 
tor. Hence the gate voltage together with the eigenener- 
gies of the DQD determine which transitions lie within 
the transport window. On the other side the value of 
the transition rates depends on the form of the eigen- 
states. In special cases a transition is forbidden since 
it violates conservation laws. A well-known example is 
the spin-blockade that occurs if the spins of the two 
states differ by more than AS = ±^ [5]. For general 
cases the effect of the structure of the wavefunctions 
enters the transition rates through the spectral weights 
[9] that determine the probability for an electron enter- 
ing a JV-electron state of the DQD to cause a transition 
to a N + 1-electron state of the DQD. 

In our calculations we assume the coupling of the 
DQD to the external reservoirs to be weak in com- 
parison with the average energy spacing in the DQD. 
Consequently the coupling to the external reser- 
voirs is treated to lowest order perturbation theory 
and the transport through the DQD is described by 
sequential-tunneling processes in and out of many- 
particle eigenstates of the isolated DQD [9]. Therefore, 
tunneling causes transitions between those eigen- 
states of the DQD which differ by one in the number 
of electrons. We assume that an electron in the up- 
per (lower) reservoir can only tunnel into the upper 
(lower) dot, which leads to different spectral weights 
for the tunneling- in and tunneling-out process [9] . As- 
suming the source reservoir to be the upper reservoir, 
the spectral weights for a transition between the i-th 
two-electron state and the j-th three-electron state 
caused by a tunneling-in (-out) process is given by 
Y2i \{N e — 2,i\di± \N e — 3, j}\ . I runs over all pos- 
sible single particle states and di+ (di-) denotes an 
annihilator of an electron in orbital I in the upper 
(lower) dot. We include the asymmetry between the 
dots caused by the source-drain voltage applied across 
the DQD by setting the asymmetry V z = ^eV s d- 
With these assumptions we calculate the stationary 



current and the differential conductance through the 
DQD in a vertical transport experiment by solving the 
rate equations containing the transition rates and the 
probabilities of occupying the many particle eigenstate 
of the DQD [9]. 

Figure 2 shows a charging diagram where the dif- 
ferential conductance G = dI/dV s d is plotted for each 
pair of transport voltage and gate voltage. The con- 
ductance shows significant values only for source-drain 
voltages Vsd > 0.35 mV. This is surprising since one 
could expect that the maxima of the conductance form 
Coulomb diamonds starting from V s d = mV and 
Vg = 16.1 mV, since at these voltages the Fermi en- 
ergies of the external reservoirs align with the energy 
needed for a GS-GS transition between two and three 
electrons. The results shown in Fig. 2 resemble the 
well-known spin blockade[5]. In fact our results can be 
explained by an isospin blockade due to the charge po- 
larization of the three electron GS. An electron tunnel- 
ing from the source reservoir (upper reservoir) to the 
upper quantum dot can hardly cause a transition from 
the unpolarized two-electron GS to the three-electron 
GS that has two electron charges in the lower dot. In 
other words, the two electron state is approximately 
an eigenstate to I z — 0, which by an electron entering 
the upper dot turns approximately to an eigenstate to 
I z = |. This state has no overlap with the three elec- 
tron GS, that is approximately an eigenstate to I z = 
— \. This explains why the spectral weight and there- 
fore the transition rate for the tunneling-in process 
nearly vanish and illuminates the relation to the spin 
blockade, which is based on the orthogonality of eigen- 
states with different spin. Rising the source-drain volt- 
age other transitions enter the transport window that 
enable transport. The strong suppression of the cur- 
rent discussed here demonstrates again the influence 
of the structure of the eigenf unctions on the transport 
characteristics [9]. 

Another prominent effect of charge localization on 
the transport properties is the appearance of broad re- 
gions of negative differential conductances in the charg- 
ing diagrams for magnetic field strengths close to the 
anticrossings. Figure 3 shows the current (upper part) 
and the differential conductance (lower part) for B = 
7.24 T. For small transport voltages, the three-electron 
GS is only weakly polarized and the GS-GS transition 
between two and three electrons is allowed. Therefore, 
Fig. 3 shows current already for small transport volt- 
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Fig. 2. Section of charging diagram showing the differential 
conductance G = 8I/dU s d for vertical transport through 
DQD. B = 7.74 T (i.e. three-electron GS is strongly po- 
larized). The conductance is strongly suppressed at low 
transport voltages, due to an isospin blockade. Units of G: 
eT/mV, where F = DOS \Tr\ 2 2zl determines the coupling 
to the external reservoirs, DOS is the density of states in the 
reservoirs, T^: tunneling matrix elements to the reservoirs. 
Parameters like in Fig. 1 and temperature T = 140 mK. 

ages in contrast to Fig. 2. But increasing the source- 
drain voltage also increases the asymmetry between 
the dots and therefore the polarization of the three- 
electron GS, which reduces the spectral weight for the 
tunneling-in process. Hence the current may decrease 
with increasing transport voltage giving rise to nega- 
tive differential conductances. 

In conclusion: We observed and explained correla- 
tion induced parity crossings in a symmetric DQD 
in a magnetic field sweep. These crossings convert to 
a magnetic-field-dependent charge-localization for an 
odd numbers of electrons inside the DQD if the sym- 
metry between the dots is broken. This localization 
strongly affects the transport characteristics and leads 
to an isospin blockade as well as to negative differential 
conductances. We thank Michael Tews for many illumi- 
nating discussions and acknowledge financial support 
by the Deutsche Forschungsgemeinschaft via SFB 508. 
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Fig. 3. Charging diagrams for transport through DQD. 
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